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ABSTRACT 



Aims. We test the isotropy of the expansion of the Universe by estimating the hemispherical anisotropy of supernova type la (SN la) 
Hubble diagrams at low redshifts (z < 0.2). 

Methods. We compare the best fit Hubble diagrams in pairs of hemispheres and search for the maximal asymmetric orientation. For 
an isotropic Universe, we expect only a small asymmetry due to noise and the presence of nearby structures. This test does not depend 
on the assumed content of the Universe, the assumed model of gravity, or the spatial curvature of the Universe. The expectation 
for possible fluctuations due to large scale structure is evaluated for the A cold dark matter (ACDM) model and is compared to the 
supernova data from the Constitution set for four diff'erent light curve fitters, thus allowing a study of the systematic eff'ects. 
Results. The expected order of magnitude of the hemispherical asymmetry of the Hubble expansion agrees with the observed one. 
The direction of the Hubble asymmetry is established at 95% confidence level (C.L.) using both, the MLCS2k2 and the SALT II light 
curve fitter. The highest expansion rate is found towards (£, b) ^ (-35°, -19°), which agrees with directions reported by other studies. 
Its amplitude is not in contradiction to expectations from the ACDM model. The measured Hubble anisotropy is AH/H ~ 0.026. With 
95% C.L. the expansion asymmetry is AH/H < 0.038. 



1. Introduction 

Isotropy of the distributions of matter and light in the Universe is 
one of the basic assumptions of modern cosmology. The small- 
ness of the temperature anisotropics in the cosmic microwave 
sky shows that this assumption is an excellent approximation. 
However, the cosmic microwave sky probes the symmetry of the 
Universe at the time of photon decoupling and it is interesting 
to test the assumed isotropy at smaller redshifts. A recent study 
of various probes at much smal ler redshift reveals fluctuations at 
the per cent level and larger (Gibelyou & H utereni201 2). These 
can partly be explained by local large scale structure, partly by 
systematics of the observations, but some of them might be at 
odds with our understanding of cosmology and call for further 
investigation. 

Complementing the assumption of isotropy with the addi- 
tional assumption of homogeneity predicts the space-time metric 
to become of the Robertson- Walker type, predicts the redshift of 
light z, and predicts the Hubble expansion of the Universe. Then 
the cosmic luminosity distance-redshift relation for comoving 
observers and sources becomes 



ddz) 



cz 



l + (l-^o)- 



-0(.3) 



(1) 



with Hq and qo denoting the Hubble and deceleration param- 
eters, respectively. Note that this prediction holds for arbitrary 
spatial curvature, any theory of gravity (as long as space-time is 



* kalus at physik dot uni-bielefeld dot de 
dschwarz at physik dot uni-bielefeld dot de 
niaj-ina dot seikel at uct dot ac dot za 

^ alexander dot wiegand at aei dot mpg dot de 



described by a single metric) and arbitrary matter content of the 
Universe. 

As suming the cosmological principle, ISeikel & Schward 
( 2009) showed that the acceleration of the cosmic expansion can 
be established in a model independent test (any curvature, any 
matter content, any dynamics of gravity) at 4.2 sigma from the 
Union data set of supernovae. One of the assumptions of this test 
is that supernovae are a fair sample of i s otrop ically distributed 
standard candles. In Seikel & Schwarj (12009 ) it was demon- 
strated that local SNe (z < 0.2) are crucial to this conclusion. 
These local SNe come from a tiny fraction of the Hubble vol- 
ume only, thus it is not obvious, if they form a fair sample and if 
they are sampled from an isotropic and homogeneous distribu- 
tion. Testing this assumption is the purpose of this work. 

A simple test for isotropy relies on comparing Hubble dia- 
grams for pairs of hemispheres. Statistically significant devia- 
tions from an isotropic Hubble diagram have been discovered by 
Schwarz & Weinhorst (2007) at redshifts z < 0.2 in three dif- 
ferent data sets. However, the direction of maximal hemispheric 
asymmetry is close to the equatorial poles, and thus a systematic 
error in the SN search, observation, analysis or data reduction 
seems to be plausible. 

A more natur al direction of anisotropy was identified by 
IColin et aP (1201 ll) . They found out that the lo w redshift regime 
(z < .06) of the data in the Union 2 catalogue (lAmanuUah et al.l 
l2010t) is barely consistent with ACDM at 2 - 3 cr because of a 
bulk flow towards the Shapley supercluster. 

Anto niou & Per ivolaropoulosI (l2010h also applied the hemi- 
sphere comparison method to the SN la data of Union 2 
and found a maximum anisotropy, which itself is statisti- 
cally not significant, but coincides with the bulk velocity flow 
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axis (iFeldman et al.ll2QlQl: iKashlinsky et al.|[2009l : iLavaux et 
12010'), co smic m icrowave background (CMB) low multipole 
moments ( iBenn ett et al. 2011; Copi et al. 2010) a nd th e quasar 
optical polarisation alignment axis (Hutse mekers et al.l 12008). 
Similar directions can be obtained by considering a dynami- 
cal dark energy fluid with an anisotropic equation of state in 
the wCDM model or in Chevallier-Polarski-Linder parametri- 
sation (Cai & Tuol 1201 2l) SNe la have also been used by 
iKoivisto & Motal (12008 alb|) to test othe r model s of anisotropic 
dark energy and by iKoivisto et al. I (l201lh to constrain 
anisotr opic curvature. Recently, ariano & PerivolaropoulosI 
(l2Q12h detected an alignment of supernova dipole asymmetries 
and a dipole of a locally varying fine structure constant, which 
can be observ ed in the Keck-hVLT quasar absorber sample. 
iJacksonI (l2012h found a hemispherical anisotropy in the mat- 
ter density Qrn at higher redshifts (0.5 < z < 3.787), where 
the smallest value of flm points close towards the CMB dipole. 
An additional tool to test the assumption of isotropy consists 
of measuring polarisations of diff'erent astrophysical sources, as 
anisotropic expansion has eff'ects in electrodynamics leading to 
observational consequences (Ciarcelluti 2012). 

In the present paper, we analyse directions of maximally 
asymmetric Hubble expansion constraining its maximum value. 
Besides, we infer from SN data typical fluctuations of the ex- 
pansion rate and compare them with expectations for the ACDM 
model. The use of S N data fitted with 4 diff'erent light-curve fit- 
ters by lHicken et al.l |2009) allows us to study systematic eff'ects. 

The paper is structured as follows: We explain the method 
of our test and state the corresponding theoretical expectation 
on the fluctuations in section [2l We continue with describing the 
data set, which we will analyse in section O Statistical analy- 
ses can be found in section |4] and possible systematics will be 
discussed in section [S] Upper and lower limits on the expansion 
asymmetry are given in section |6l We conclude in section [7] 



2. Hemispherical asymmetries 

Our aim is to test whether the Hubble diagram is isotropic. 
In orde r to do so, we apply a simple test, that first has been 
used by lEriksen et al.l (|2004|) to study asymmetr i es in the CMB 
anisotropy field, and that lSchwarz & WeinhorstI (12007) adopted 
to find anisotropics in SN data. We test the anisotropic cosmic 
expansion by splitting the sky into hemispheres. Having calcu- 
lated the Hubble diagrams on both hemispheres separately, one 
can change the directions of the poles and recalculate the Hubble 
diagrams on the new hemispheres. 

We fit the Hubble diagrams by applying the method of least 
squares to the distribution of distance moduli of the supemovae, 
using as a model the relation of the distance modulus to the red- 
shift given by ([T]) and the definition 



jd = 5 logio (dL/Mpc) -h 25 . 



(2) 



ISchwarz & Weinhorsi (l2007h estimated, that the error from trun- 
cating Hubble's law as in ([T]) is below 10% if one fits only at 
redshifts z < 0.2. Despite this restriction, this test has the major 
advantage, that it is model-independent, because ([B follows di- 
rectly from the cosmological principle and does not contain any 
terms depending on the curvature, theory of gravity or content 
of the Universe. 

Admittedly, as to go from SN magnitudes m = yu -h At to 
the luminosity distance di, one has to know the nuisance pa- 
rameter M = M - 5 log 10 (^o) + 25, which also depends on 
the absolute magnitude M and hence on the physics of SNe 



(iPerlmutter et al .1119991) . We assume that Type la SNe are stan- 
dardisable, thus we trea t At as a universal number and use the 
distance moduli which iHicken et al. ( 2009) calculated by us- 
ing four diff'erent light-curve fitters with calibrations Hq = 65 
km/s/Mpc and Mb = -19.46 for SALT, Mb = -19.44 for SALT 
II and My = -19.504 for MLCS2k2. 

When fitting the Hubble diagrams, we allow Hq to vary from 
hemisphere to hemisphere, but keep the deceleration parameter 
at its best-fit ACDM model value. The SN Hubble diagrams con- 
sidered here do not constrain qo well enough, to justify the in- 
clusion of this parameter in the fit. 

In a perfectly homogeneous and isotropic Universe we 
would expect any nontrivial result in the hemispherical asym- 
metry test to be due to noise and systematics. However, 
the large scale structure of the Universe is expected to give 
rise to a small deviation from isotropy. We can estimate 
the expected amount of fluctuations of the Hubble rate be- 
twee n the hemispheres by estimating their cosmic variance 
(Tur ner etal.1 1 19921; Ishilll997l; IWang et al.l [19981: iLi & Schwarzl 
2008?Wiegan d & Schwar3l2012 h. To measure the diff'erence be- 
tween the two hemispheres we use the quantity • In the for- 
malism o f IWiegand & Schwarzl (l2012h the Hubble rate on each 
of the hemispheres is given by the volume average over the lo- 
cal expansion rates. In the linear regime, these are related to the 
density contrast and therefore the average Hubble rate on a hemi- 
sphere is 



S/N 



S/N 



1 



:fvo {^o}s/N 



(3) 



where the growth rate factor today f^o^ ~ 0.5 in the standard 
ACDM model. To first order the anisotropy therefore is 



(4) 



So the ensemble expectation value of the anisotropy is zero at 
leading order. However, we are interested in the typical fluctua- 
tions between difl'erent locations in the Universe. Using (|4]) we 
find 



cta 



Hm — Hs 



HS 



FS 



where 



is the matter variance on a hemisphere and 
oris •= «^o)d) = f Po (k) Wl (k) d'k . 

J]R3 



(5) 



(6) 



(7) 



is the matter variance on the full sphere. Pq (k) is the linear mat- 
ter power spectrum and W^o (k) is the Fourier transform of the 
radial window function W^) (r) of the domain D. This window 
function is often taken to be of top-hat- or Gaussian-form. To 
arrive at more precise predictions, however, we shall use in Sec. 
[3] a shape that describes the actual SN distribution more accu- 
rately. 

The explicit form of cr^^ is more complicated than the one 
for the full sphere, because a window function with the geometry 
of a hemisphere is no longer spherically symmetric. Therefore, 
the window function has also an angular part in addition to the 
radial part (r). We calculate this angular part using a decom- 
position into spherical harmonics (see appendix A of iWiegandl 
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(l2012h for the explicit form). Apart from this angular window 
function, which is fixed by the requirement that we consider a 
hemisphere, cr^^ also uniquely dep ends on Pq (k) SindW^ (r). 
Thus, using the fitting formula of lEisenstein & Hul (Il998 ) for 
the matter power spectrum allows us to estimate the possible 
anisotropics in a ACDM universe in Sec. [3] 

3. Supernova data and results 

For our analysis, we use the Constitution set (iHicken et al.ll2QQ9b 
supplemented by the positions of the SNe procured from the list 
of SNe Q provided by the lAU Central Bureau for Astronomical 
Telegrams (CBAT). The Constitution set contains a large num- 
ber of nearby SNe from all the sky (except for the zone 
of avoidance), whose distance moduli have been obtained by 
four diff'erent light-curve fitters: Multicolor Light Curve Shapes 
2k2 (MLCS2k2) (1.7), MLCS2k2(3.1), Spectral Adaptive Light 
curve Template (SALT) and SALT IL Thus, the Constitution set 
consists of four s amples. 

The SALT (Guv et al.|[2QQ5 ) fitter describes Hght-curves by 
the peak apparent 5-band magnitude m^^^, the stretch factor s 
and the color parameter c. SALT does not distinguish between 
the intrinsic color of a SN and the reddening which is caused by 
dust extinction. Both eff'ects are summarised in c, thus assuming 
the same color-magnitude relation for the combination of color 
variation and hos t reddening at low and at high redshift. SALT II 
(iGuy et al .1120071) is similar to SALT, but uses a diff'erent spectral 
template, s has been replaced by a new stretch param eter xi . 

In contrast to the SALT fitters, MLCS2k2 (Jha etalJl2007|) 
includes the intrinsic color in the shape parameter A assuming 
a broader-bluer and narrower-redder relation. The host-galaxy 
extinction parameter Ay is then determined separately. This re- 
quires assumptions about the dust properties specified by the 
reddening parameter Ry. In the Milky Way, this parameter has 
been measured to be = 3 . 1 . We will r efer to the MLCS2k 2 fit- 
ter using this value as MLCS2k2(3.1). Hicken et all (l2009l) have 
determined the reddening parameter, which minimises the scat- 
ter in the Hubble diagram, to he Ry = 1.7. In the following, the 
fitter with this value will be denoted as MLCS2k2(1.7). 

While there is a significant overlap of the SNe contained in 
the samples, not all observed SNe can be found in a specific sam- 
ple. SNe are rejected due to certain criteria, e.g. if the light-curve 
fit is not satisfactory or if the light extinction is too high. The 
diff'erent parameters and rejection criteria of the fitters thus lead 
to a diff'ering number of SNe in the samples. Fig. [T] shows the 
distribution of SNe on the sky for all four samples. In Table [T] 
the number of SNe with z < 0.2, their mean and median red- 
shifts and their mean cr^ are listed for all four samples. The 
SALT sample contains a significantly smaller number of SNe 
than the other three sa mples as the stricter rejection criteria from 
iKowalski et al.l ([2008) have been adopted. 

In our data analysis, we set the deceleration parameter qo = 
-Qa -h = -0.601 according to the seven-year Wilkinson 
Microw ave Anisotropy Pro be (WMAP) observations fitted to 
ACDM (iLarson et al.ll20TTh . 

The eff'ect of our choice of qo is not too large: Repeating our 
test described below for SALT II and fixing qo at diff'erent val- 
ues between -1.5 and 0.4, yields values of (Hn - ^5)max which 
vary almost linearly between 3.20 km/s/Mpc for = -1-5 and 
(Hn - ^5)max = 2.55 km/s/Mpc for qo = 0.4. 

As said in section[Tl we carry out a;^^-fit to adjust the Hubble 
parameters Hj^ii^b) on the hemisphere identified by its pole 

^ http : //www . cbat . eps . harvard . edu/lists/Supernovae . htir 




Fig. 1. Celestial distribution of SNe from 'Hicke n et al] (l2009b 
having redshifts z < 0.2 in galactic coordinates. Black dots rep- 
resent SNe, that are in all four samples, while green ones can- 
not be found in the SALT sample, blue ones are only in the 
MLCS2k2 samples. The MLCS2k2 (Ry = 3.1) sample is the 
only subset, which contains red SNe. We plot in Mollweide pro- 
jection. 



fitter 


# SN with 


mean 


median 


mean cr^ 




z<0.2 


redshift 


redshift 


in mag 


MLCS2k2(L7) 


199 


0.030 


0.025 


0.191 


MLCS2k2(3.1) 


203 


0.030 


0.025 


0.197 


SALT 


115 


0.037 


0.028 


0.187 


SALT II 


183 


0.031 


0.026 


0.218 



Table 1. Compilation of some characteristics of the samples 
used throughout this article. 
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Fig. 2. Colour scaling of all maps drawn for this article, - 
Hs denotes the difference in the expansion rates for a pair of 
hemispheres, is the corresponding number asymmetry 

of SNe, is the asymmetry in dust antenna temperature, 

stands for the asymmetry in the shape parameter of the 
MLCS2k2 fitter, and is the colour parameter asymmetry 
(SALT/SALT II). 

(£, b) and Hg if, b) on the complementary hemisphere such that 
the distance moduli modelled by ^ a nd related to redshif t data 
by ([B fits the distance moduli given in lHicken et al.l (|2009^ best. 
We divide the sky into a 1° x 1° grid. Each grid point can be 
regarded as the pole (^, b) of a hemisphere. We determine the 
Hubble parameter Hj^^i, b) and subtract from it the Hubble pa- 
rameter Hs if, b) of the opposite side. The resulting deviations in 
Ho are plotted for each light-curve fitter in Fig. [3] 

The asymmetry maps for the SALT fitter exhibits several 
small local extrema in Hf^ - Hg , which look like noise, while 
the maps obtained with the MLCS2k2 fitters and SALT II are 
similar to each other: The white areas, which denote small devi- 
ations in the Hubble parameter, have similar shapes, and the di- 
rections of maximum asymmetry are also similar, as can be veri- 
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Fig. 3. Hemispherical asymmetry - Hs in the Hubble rate for each SN la fitter. The upper left plot is for MLCS2k2 (1.7), 
the upper right one is for MLCS2k2 (3.1), below are SALT (left) and S ALT II (right). The deceleration parameter is fixed at the 
value qo = -0.601, which one obtains from WMAP ( Larson et al.ll201 ll ). When the north pole lies in a region with a bluish colour, 
- Hs is negative, the red spots denote the directions of positive large asymmetry - Hs . The Hubble parameter has no 
anisotropy pointing towards white areas. All plots showing the asymmetry in the Hubble rate have the same scaling. The scale 
is given in Fig. O At the reddest and bluest spots, \Hn - Hs\ = 3.4. We plot directions in galactic coordinates and Mollweide 
projection. The black line is the equator of the equatorial coordinate system. 



fied in Tabled Table [2] also shows, that the maxima of ^^"^^ for 
the MLCS2k2 fitters and SALT II are well above 2% and hence 
larger than the fluctuations cr^ due to cosmic variance (cf . end of 
this section). MLCS2k2(1.7) and SALT II yield almost the same 
directions of maximum asymmetry. The maximally asymmetric 
direction in the MLCS2k2(3. 1) is located in the same quadrant as 
the previous ones, whilst the direction found in the SALT data 
is far from the others. One finds also a much smaller value of 
§^ = 1.54% using the SALT data. 

For a correct interpretation of these anisotropics, we need to 
compare them to what we expect by cosmic variance. In a sta- 
tistically homogeneous and isotropic ACDM model, the average 
hemispherical asymmetry is zero. A particular realisation, how- 
ever, will have an asymmetry. The broadness of the distribution 
of these asymmetries around zero, is the cosmic variance for the 
hemispherical asymmetry. We use © to estimate it for the four 
samples. To arrive at reliable predictions for the expected fluctu- 
ations in the anisotropy, we have to account for the incomplete- 
ness of the sampling of the hemispheres by the supemovae. To 
this end, we use a radial window function that is adapted to the 
shape of the redshift distribution of the SNe. We show this red- 
shift distribution in Fig. [4] together with the results of a fit with 
the ansatz Ar^/ (l -\- 5r^^. We use this shape as the radial win- 
dow function (r) oc A/ (l -\- Br"^^ in the evaluation of Q and 

Plugging the resulting values for cr^^ and cr^^ into ([5]), we 
therefore arrive at a value of about 0.7% as shown in Tabled 
Note that this is the expectation value for a randomly fixed direc- 
tion, not for the direction of maximum asymmetry. We checked 



fitter 




b 


(H^ - Hs )tnax 


Hn-Hs 
Hn+Hs 


MLCS2k2(1.7) 


-52° 


-19° 


3.3 km/s/Mpc 


2.49% 


MLCS2k2(3.1) 


-18° 


-20° 


3.4 km/s/Mpc 


2.59% 


SALT 


-154° 


32° 


2.0 km/s/Mpc 


1.54% 


SALT II 


-52° 


-18° 


3.0 km/s/Mpc 


2.28% 



Table 2. Galactic longitudes ^ and latitudes b of the directions 
of maximum asymmetry (H^ - Hs )max in the Hubble rate ifo if 
the deceleration parameter qo is fixed. 



if the shape of the radial window function is important by using 
diff'erent forms, including a simple interpolation of the distribu- 
tion of SN. The values varied between 0.6% and 1.0% indicating 
that, in any case, the cosmic variance is well below 2%. 

On the other hand. Table [T] tells us, that the mean measure- 
ment uncertainty cr^ of the distance modulus is approximately 
0.2 mag and that the number of SNe with z < 0.2 is around 
200 for the MLCS2k2 data. Thus, we calculate the Hubble pa- 
rameters with roughly 100 objects per hemisphere, reducing the 
overall uncertainty in jd to an approximate value of 0.02, which 
relates to an uncertainty in the Hubble parameter of the order of 
0.01. 

Hence, the sensitivity of this test applied to the Constitution 
set is at the same order as the expected fluctuations of ACDM, 
which we therefore cannot see. However, more extreme models 
would stick out. 
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Fig. 4. The radial distribution of the supemovae in the MLCS2k2 (1.7) (top left), MLCS2k2 (3.1) (top right), SALT (bottom left) 
and SALT II (bottom right) data set. The red dashed lines illustrate the radial window functions which we use in our calculations 
for the expected fluctuations. 



4. Statistical analysis 

Having found asymmetries, we want to focus on the question of 
the likeliness of the maximal asymmetry observed on the sky. 
Therefore, we performed the following analyses: 



4.1. Scrambled data 

In our first approach we do not relate a SN from the Constitution 
set to its original position from the CBAT list, but relocate it 
at the position of another random SN, that is also part of the 
same sample of SNe. We repeat the procedures described in 
the previous section for 16000 realisations, which were gen- 
erated by scrambling the positions. Thereafter, we again max- 
imise the asymmetries in - Hg, but as to save computing 
time, we do not calculate the asymmetry in every grid point but 
first evaluate the hemispherical Hubble asymmetry at 10 diff'er- 
ent random positions and then use the maximum of these as 
a start po sition (^,^7)0 for th e simulated annealing method de- 
scribed in iHag gstromi fe002') . In this method a Metropolis chain 
is simulated by determining the Hubble asymmetry on a ran- 
domly picked neighbouring point (£,b)n+i = (^„ ± 1°,Z7„) or 
(^, b)n+i = i^n^i^n ± 1°)- The new position (^, b)n+i is either ac- 
cepted with a probability 



mm^ 



[exp( 



1 H^-Hs 



,1 



(8) 



or rejected with 1 - /?acc, i-e. (£,b)n+i = (^,b)n. The algorithm 
converges to the position lim„^oo(^, /?)« = (^, /?)max at which 
the hemispher ical Hubble asymmetr y is maximised, as long as 
Tn oc 1^ (Gem an & Gemanlll984l) . The algorithm is trained 
to find the same Hubble asymmetry values as in Table |2] with 
the same accuracy. Out of 10 runs, we identified the maximum 
number of steps needed to obtain these values and we chose 
a slightly higher number of steps, after which we truncate the 
Metropolis chain and save the highest expansion asymmetry ap- 
pearing therein. 

Assuming a perfectly isotropic Hubble expansion, any devia- 
tion from hemispherical symmetry would be due to selection ef- 
fects. Thus scrambling the original data yields an estimate of the 
distribution one would expect for an isotropic universe, which is 
measured with the true selection function. Moreover, no model 
has to be adopted for this test. 

The results are visualised in the histograms of Fig. [5] There, 
one can see again a distinction between the results for the SALT 
and all the other data samples: the maximum asymmetry in Hq 
for the SALT data set is a fair realisation of this Monte Carlo 
procedure, whilst between 90 and 95 % of the maxima of all 
Monte Carlo realisations are smaller than the asymmetry in the 
Hubble diagrams for the other light-curve fitters. 

4.2. Distance modulus Monte Carlo 

As we filled in the previous analysis observed positions with 
observed SNe, the simulated data can also be aff'ected by an 
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Fig. 5. 4 Histograms showing the maximal hemispherical asymmetry in the Hubble rate Hq after scrambling 16000 times the 
positions of SNe in the MLCS2k2 (1.7) (top left), MLCS2k2 (3.1) (top right), SALT (bottom left) and SALT II (bottom right) 
data sets. Corresponding to Fig. [3] the deceleration parameter qo = -0.601 is fixed. The red line denotes the value one gets by 
analysing the original data from lHicken et al.l (l2009l) . the dashed lines denote the quantiles. 



anisotropic distribution of observed SNe. Another way to check 
whether the asymmetries are significant is to repeat the follow- 
ing procedure: 

1 . We generate a data set, in which we keep the positions and 
redshifts of each SN, but the distance modulus jd is cal- 
culated from its redshift z according to (O, where we as- 
sume the WMAP A CDM values = -0.601 and Hq = 
71.0 km/s/Mpc from iLarson et aP (1201 ih . 

2. We simulate "observations" with uncertainties of measure- 
ment as a Gaussian distribution N(jd, cr) , where the expecta- 
tion value jd equals the previously calculated distance mod- 
ulus jd and the standard deviation cr equals the measurement 
uncertainty cr^. of the distance modulus of the ith SN in the 
original data set. 

3. We maximise again the anisotropy - Hs in the simu- 
lated data set by dint of the simulated annealing method 
(lHaggstromll2002h . 

As the Hubble parameter is normally distributed, and thence 
the diff'erence - Hs as well, we expect the maximum asym- 
metry in Hq to obey a generalised extrenie valu e di stributioifl 
due to the Theorem of Ipisher & Tippetj (1 19281) and iGnedenkol 

^ Let Xi,X2,X2, ... be independent and identically distributed (iid) 
random variables. Let F be the underlying cumulative distribution 
function. Then the probability, that the maximum of the iid ran- 
dom variables is smaller than x, is P(max (Xi, ...,X„) < x) = F^{x). 
Suppose there exists a sequence of constants Un > 0, and bn e 

has a non-degenerate limit distribution 



such that 



an 



fitter 




cr 




l-CDF(f^) 

\H^+Hs /max 


MLCS2k2(1.7) 


2.11 


0.47 


-0.10 


5.1% 


MLCS2k2(3.1) 


2.14 


0.48 


-0.12 


4.4% 


SALT 


2.40 


0.56 


-0.11 


86.2% 


SALT II 


2.47 


0.55 


-0.11 


30.3% 



Table 3. Maximum likelihood parameter estimates (yu, cr, ^) ob- 
tained from the distance modulus Monte-Carlo with diff'erent 
light-curve fitters (cf. figO. The last column contains the proba- 
bility of measuring higher Hubble asymmetries than the original 
ones according to the fitted distribution. 



(1 19431) . The probability density function of which with location 
parameter jd, scale parameter cr and shape parameter ^ 9^ is 
given by 



f(x;/d,cr,^) = - 
cr 



1 -h 



-i/f-i 



1/^ 



(9) 



The maximum likelihood parameter estimates can be found in 
Table [3] For an infinite number of realisations, we would expect 
^ = 0, also known as Gumbel distribution. As shown in Table [3j 
1^1 « 1. 

The resulting histograms (Fig. (6]) show statistical signifi- 
cance for the anisotropics found in the MLCS2k2 data, which 



lim^^^oo (UnX + bn) = G( x). Such a distribution G is called an extreme 
value distribution (see e.g. lde Haan & Ferreiral(l2006h ). 
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Fig. 6. Histograms of the maximum asymmetry in ifo found in 16000 different realisations of the distance modulus Monte Carlo, 
applied on the MLCS2k2 1.7, MLCS2k2 3.1, SALT and SALT II data sample, and the respective probability density functions of a 
generalised extreme value distribution. The distance moduli in the simulated data are normally distributed around their theoretically 
calculated values with the original measurement errors cr^. as standard deviation. 



are higher than 95% of the anisotropics in the simulated data, but 
they show no significance for the SALT and SALT II asymme- 
tries. The distribution for MLCS2k2 is narrower than the others, 
because the SALT sample contains less SNe and because SALT 
II yields larger errors for the distance moduli. 

As to study how improving measurements of aff'ects the 
usefulness of our hemispheric asymmetry test, we repeat 5000 
distance modulus Monte Carlos, but we replace cr^ by (T^/2. We 
can so show, that if the measurement uncertainty of the distance 
moduli of every single supernova can be halved, the Hubble 
asymmetries are expected to be much smaller. Fitting again 
a generalised extreme value distribution, the maximum likeli- 
hood location parameters would amount 1.09 km/s/Mpc and 
1.24 km/s/Mpc for MLCS2k2 (3.1) and SALT II respectively. 
The observed values would lie at the 1 - 10"^ and 1 - 4 x 10"^ 
quantiles for the MLCS2k2 (3.1) and SALT II Hght-curve fitters. 
We hope, this will encourage observers to further reduce mea- 
surement errors, as it also improves this test considerably. 

4.3. Comparison with cosmic variance 

Finally, we want to determine the fluctuations of the hemispher- 
ical asymmetry around its average, for a random orientation of 
the hemispheres. To calculate these fluctuations, we pick 10,000 
random positions (£,b) in the sky. Each position so determines 
the north pole of a random sphere. We split this sphere into 
its northern and southern hemisphere and calculate from the 



Constitution set the hemispherical Hubble asymmetry ^^^^^ of 
this randomly oriented sphere. The variance of the Hubble asym- 
metry for the 10,000 spheres is then the desired fluctuation. 

To sample the random north pole positions (£,b), we employ 
a restriction: If we chose them from the whole sky [£ uniformly 
distributed, pdf(Z7) = | cos(Z7)], for each position where we take 
(£,b) to be the north pole of our random sphere, we would also 
have a random sphere in the sample where (£,b) is the south pole, 
because we chose the position (tt -h £,-b) as a north pole. So in 
order for the average to be the typical asymmetry for random 
orientation, we have to restrict the choice of random north pole 
positions (£,b) to one hemisphere of the sky. 

In order not to miss a direction of high anisotropy, we do 
three runs for orthogonal directions in the sky that we take as 
the poles of the hemispheres on which we sample our random 
positions (£,b). In the first run, we generate only positions (£,b) 
on the northern galactic hemisphere, in the second run positions 
are restricted to the western hemisphere and in the third repeti- 
tion the generated (£,b) are located on the hemisphere towards 
the galactic centre. 

Table |4] compares the empirical root mean square fluctua- 
tions (Temp of the hemispherical Hubble asymmetry determined 
by this prescription with the fluctuations in the Hubble asymme- 
try due to cosmic variance cta . Under the assumption that each 
configuration of SN in the sphere obtained by the random choice 
of its north pole position represents a realisation of the stochastic 
process that generates the SN distribution, (Temp should give an 
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fitter 



hemisphere 



( Hn-Hs \ 
'\Hn+Hs j 



O-A 



SALT2 lAHol 

l.Or 



ATT /^COI^O /I '7\ 

MLCoZKZ (i. /) 




O.qVc 


f\ noi 




MTT 


1 f\f7/ 


c\ no/ 
U. r/c 




Wrl 


1 no/ 
L.UVc 


rv no/ 


MLL^b ZkZ (3.1) 




rv no/ 


0.7% 




NH 


0.8% 


0.7% 




WH 


0.8% 


0.7% 


SALT 


GC 


0.6% 


0.7% 




NH 


0.5% 


0.7% 




WH 


0.6% 


0.7% 


SALT II 


GC 


0.8% 


0.6% 




NH 


0.9% 


0.6% 




WH 


0.7% 


0.6% 



Table 4. Empirical standard deviation (Temp of the distribution of 
evaluated on random positions on the northern (NH) and 
western (WH) hemisphere as well as on the hemisphere towards 
the galactic centre (GC) compared with expected fluctuations a a 
due to cosmic variance. 



estimate of a a - Table |4] shows, that their values indeed are quite 
comparable. Only the empirical values obtained with MLCS2k2 
(1.7) and SALT II are slightly higher than expected from this 
reasoning. 



5. Systematics 

5.1. Equatorial hemispheres 

In this section we are concerned with the question how un- 
likely is the asymmetry observed towards a fixed direction on 
the sky. Therefore, we carry out again the distance modulus 
Monte Carlo, but now we do not maximise the anisotropy, but 
we calculate the diflTerences in the Hubble parameters of the 
equatorial hemispheres. The results all look like those obtained 
from SALT II (Fig. (T]), where the measured value Hj^ - Hs = 
1.15 ^ = = 0.88% lies within the 90% confidence 

interval of the resulting distribution. The corresponding values 
for MLCS2k2 1.7 and 3.1, as wefl as for SALT are 1.00 

' s-Mpc ' 

1.25 ^ and -0.45 ^ or accordingly 0.76%, 0.96% and 
-0.34%. These values are all smaller than the cosmic variance. 
Thence and because the values lie within their corresponding 
90% confidence intervals, this test does not support that the ob- 
servation strongly depend on Earth related effects. 



5.2. Number asymmetry 

As already mentioned, the asymmetries in can be due to in- 
complete SNe la observations in some directions. Therefore, we 
plotted the diflference between the number of SNe on both hemi- 
spheres similarly to Fig. [3] These plots (Fig. [8]) show anisotropics 



in 



Nn-Ns 
Nn+Ns 



but they do not agree with any asymmetry map in Fig. 
[3l The asymmetry in the SN distribution is minimal near the 
Earth's equator and maximal approximately at its poles. Also 
the white region, in which the number asymmetry is zero, fol- 
lows the equatorial equator. This is not the case in Fig. [3] 

Although the maximal asymmetries in the number of SNe 
is around 50% in all data sets, the number asymmetries for 
the hemispheres that give the maximal expansion anisotropics 
(cf. Table O are much smaller for MLCS2k2 and SALT II 
(MLCS2k2 1.7: -4.8%, MLCS2k2 3.1: -20.2%, SALT II: - 
10.6%). Only SALT manifests a comparably high number asym- 




Fig. 7. Histogram with the resulting diflTerences between the 
Hubble rate on the northern equatorial hemisphere and the 
southern equatorial hemisphere Hs obtained from 16000 runs 
of the distance modulus Monte Carlo. The red line indicates the 



measured value Hm - Hs = 1.15 



km 



obtained from the SALT 



iN-ns --s-Mpc 
II data set. We also plotted a normal distribution N(pi, cr) with 
yu = 0.00 and cr = 1.07. 




Fig. 9. Hemispherical asymmetry in the antenna teniperat ure of 
the foreground dust extracted from iSchlegel et al.l (Il998l) . The 
maximum asymmetry is 0.14 mK. 



metry between the hemispheres of maximum expansion asym- 
metry, that is to say 46.9%, even though it is the most isotropic 
subsample in terms of Hubble asymmetry. 

For that reason we cannot think of the incompleteness of the 
SN data being the main explanation of asymmetries found in the 
Hubble expansion. 



5.3. Quality of fit 

Another possible systematic eflfect could be, that the quality of 
our fits is diflTerent in diflTerent directions. We have verified, that 
the asymmetries in ^ are also not aligned with the asymmetries 
inifo. 



5.4. Dust asymmetry 

We assume, that the Universe is isotropic on large scales. The 
foreground, like dust within our galaxy, is not isotropic and 
therefore a possible explanation of apparent anisotropic expan- 
sion. 

For the comparison with galactic dust, we brought the Schlegel 
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Fig. 8. Hemispherical asymmetry ^^^^^ in the number Ni of SNe 
same way as Fig.O The maximum value, which defines the scale, 

dust map (ISchlegel et al. into the same form as our other 

maps. The result can be seen in Fig.O This map has no similar- 
ities with Fig. [3l thus, we can also exclude dust as a reason for 
apparent anisotropic cosmic expansion. 

5.5. Shape asymmetry 

The MLCS2k2 fight -curve fitter uses the light-curve shape to 
r emove t he influ ence of dust from the distance measurement 
(iRiess etal.lll998l) . As iHicken et al.l (120091) provides the shape 
parameter of each SN, we also plot the weighted difl'erence 
^^"^^ between the average shape parameters A of each hemi- 
sphere (see Fig. [T0|. These plots also do not reproduce the 
shapes of the asymmetry map of Hq (Fig.O. 

5.6. Colour asymmetry 

SALT and SALT II u se a colour parameter to get rid of dust 
(iGuy et al.ll2005ll2007 '). We tested the asymmetries in the colour 
parameter in the same way as those in the shape parameter. The 
resulting maps are shown in Fig. [TT] The SALT data set shows 
directions in which the colour parameter is obviously asymmet- 
ric, as against irregular asymmetries of the colour parameter in 
the SALT II data. 

6. Limits on expansion asymmetry 

As to find limits on the Hubble expansion asymmetry of the lo- 
cal Universe, we draw realisations from the SN magnitudes and 
redshifts given the reported measurements with errors. We then 
compute the asymmetry of the original direction of maximum 
asymmetry given in Table [2l After 5400 realisations, we enu- 
merate the 95% quantile Hgs, which can be regarded as an upper 
limit on the expansion asymmetry, as well as the 5% quantile H5 . 
In an isotropic universe, all measurements of the Hubble asym- 




on the Hemisphere labelled with /. The plots are adjusted in the 
is 0.52. 



fitter 


5% quantile 


95% quantile 


MLCS2k2(1.7) 


0.8% 


3.4% 


MLCS2k2(3.1) 


1.1% 


3.8% 


SALT 


0.0% 


3.0% 


SALT II 


0.7% 


3.7% 



Table 5. Lower and upper limits on the hemispherical expansion 
asymmetry 



metry should agree with AH/H = at every position. Hence, a 
positive H5 cannot be explained by statistical fluctuations alone. 

H5 and Hgs can be found for every light-curve fitter in 
Table \5\ SALT is consistent with an isotropic local Universe, 
whilst the other light-curve fitters provide evidence for real 
anisotropics. MLCS2k2 (3.1) features an ^5, which is even com- 
parably high as our expectations for typical fluctuations due to 
large scale structure. 

In Fig. [T2I we plot the contours of the 5% quantiles of the 
MLCS2k2(3.1) and SALT II as error contours of the direction 
of maximum Hubble asymmetry. Both data sets agree with the 
expansion being fa stest into the direc tion of the WMAP cold 
spot I reported by Bennett et al.l (1201 ik . 

The Hgs values provide upper limits on the anisotropy of the 
Hubble expansion. The least stringent upper limit is found from 
MLCS2k2 (3.1), namely < 0.038 at 95% C.L. 

7. Conclusion 

We analyse hemispherical asymmetries AH/H = (H^ - 
Hs ) /(Hn -\-Hs) in the Hubble expansion rate at red shifts z < 0.2 
in SN data of the Constitution set published by iHicken et~an 
(2009). The observed hemispherical asymmetry of the Hubble 
expansion is in agreement with the order of magnitude that can 
be expected by cosmic variance in a ACDM universe, but is in- 
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Fig. 11. Asymmetry 



of the colour parameter. The maximum value amounts to 0.64. The plot on the left hand side shows the 



colour asymmetry in the SALT data, the right one in SALT II. 



■ JVlIjCS2k2 (3.1) 

♦ SALTD 

□ AntfHiiCKi & Ptm ropou los (201 0)( z<1.4) 

A (^i&Tuo(2012Xz<0.2) 

o ai&Tiio(2012)(z<1.4) 

T Qilmetiil.(2011)(z<0.2) 

A IVb lia no & Pe riv ola ropmilos ( 2012 ) 

♦ Schwa rz & Weinhoist ( 2007 ) data set A 
« Schwa rz & Weinhorst (2007 ) data set B 

♦ Sth wa re & Weinbt»rst (2007 ) da ta set D 




Du^i-tions of GVBa nisotropy : 

CVIBdJpole momeat (Linewea veret aL 1996) 
■k coldspotldetectedtnWMAPdata (Bennett etaL2011) 
o odd spot lldeteded in WM\Pda ta ( Chiz et aL 2007 ) 



Other Directions: 

V Sha pley supeiduster (Raychaudhu ry 1 
S Sou th Equ atoiiai Pole 



Fig. 12. 95 % confidence level contours of the SALT II (green lines) and MLCS2k2(3.1) (red lines) maximum asymmetry directions 
(MLCS2k2(3.1) ■, SALT II • ) as well as 90% confidence level contours (filled green an d filled re d respectively, filled dark 
green where they coincide) compared with directions of anisotro pies in SN data fron i Antoniou & P er ivolaropoul os (20 li)) (z < 
1.4) □, Cai&Tuo (2012)fe < 0.2 A, z < 1.4 o), Cohn et al.l OOTll) (z < 0.2 t). iMariano & PerivolaropoulosI (l2012l) A and 
Schwarz & W einhorst (2007) (data set A B D 4) , and with the direct ion of the CMB dip ole moment 9 (*Linewe aver et al.l 
1996h. the S haplev supercluster v ( Ravchaudhurvll 1 989h . the cold spots I * (iBennett et al.ll201l[) and II detected in WMAP data O 
(ICruz et al.ll2007) . and the South Equatorial Pole S . 
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consistent with the assumption of perfect isotropy and homo- 
geneity. 

However, the Hubble asymmetry in the direction of maximal 
asymmetry as identified by the MLCS2k2 light curve fitter, is at 
95% C.L. larger than in a random realisation of distance mod- 
ulus observations. All light-curve fitters, bar SALT, establish 
their respective asymmetric directions at 90% C.L. in terms of 
data scrambling. Fig. [T2l shows their ag reement with directions 
reported by Cai & Tuo (2012); Mar iano & PerivolaropoulosI 
(12012) and .Schwarz & Weinho rst (2007), which are in the vicin- 
ity to t he WMAP cold spot I iBennett et^ l201ll) . llnoue & Sil5 
(l2007h explained observed large-angle CMB anomalies with 
a pair of local dust-filled voids. These would cause both a 
cold spot in CMB data and fluctuations in the locally mea- 
sured Hubble constant as large as 2-4%. Our measured Hubble 
anisotropy of AH/H ~ 0.026 agrees with this finding. Despite 
the fact that it is significant, we so conclude that the anisotropy 
does not contradict global isotropy, because firstly it can be ex- 
plained by local structure, secondly we only analysed local SNe 
and thirdly it matches our predictions of typical fluctuations in a 
ACDM background. 

Admittedly, our test is not sensitive enough to see if the 1 % 
asymmetry detected for random hemispheres is due to ACDM- 
fluctuations of the local structure. So it should be repeated, when 
a full sky survey comprising approximately 2000 SNe becomes 
available, or when the errors of individual SNe could be reduced 
by about a factor of 2. But nevertheless we can constrain the 
expansion asymmetry of the local Universe to be less than 3.8% 
at 95% C.L. using most conservatively the results obtained with 
the MLCS2k2 (3.1) Hght-curve fitter. 
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